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Abstract. We present BeppoSAX results on the nearby 
starburst galaxy NGC 253. Although extended, a large 
fraction of the X-ray emission comes from the nuclear re- 
gion. Preliminary analysis of the LECS/MECS/PDS ~ 
0.2-60 keV data from the central 4' region indicates that 
the continuum is well fitted by two thermal models: a 
"soft" component with kT ~ 0.9 keV, and a "hard" com- 
ponent with kT ~ 6 keV absorbed by a column density of 
~1.2 x 10 22 cm~ 2 . For the first time in this object, the Fe 
K line at 6.7 keV is detected, with an equivalent width of 
^300 eV. This detection, together with the shape of the 
2-60 keV continuum, implies that most of the hard X- 
ray emission is thermal in origin, and constrains the iron 
abundances of this component to be ~0.25 of solar. Other 
lines clearly detected are Si, S and Fe L/Ne, in agreement 
with previous ASCA results. We discuss our results in the 
context of the starburst-driven galactic superwind model. 



Key words: Galaxies: individual: NGC 253 - Galaxies: 
spiral - Galaxies: starburst 



1. Introduction 

In the local universe, starburst galaxies (SBGs) are galax- 
ies with current (for typically 10 8 yr: Rieke et al. 1980) 
star formation rates significantly higher than the galactic 
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average. In the early evolution of galaxies (zft 1), star for- 
mation was certainly much higher than at present (see e.g. 
Madau et al. 1996; Hogg et al. 1998): a starburst phase 
was, therefore, normal in the early evolution of galaxies. 
Part of the interest in local SBGs stems then from the 
possibility that they resemble normal galaxies when they 
were at zlt 1. High-energy phenomena highlight the SB 
activity (Bookbinder et al. 1980): because of obscuration 
of the optical emission (by star-heated IS dust), X-rays 
allow a penetrating view of SBGs. 

Aiming at studying the X-ray properties of local SBGs, 
we have undertaken a program to study two such sources, 
NGC 253 and M82, with BeppoSAX: in this Letter the 
broad-band spectrum of the nuclear region of NGC 253 is 
reported, focusing on the discovery of the Fe K line at 6.7 
keV and on the features of the continuum in the 2-10 keV 
band arising from the nuclear region. 

Assumed to be a prototypical SBG, NGC 253 has been 
extensively observed in X-rays. At low energies, analysis 
of Einstein IPC data led Fabbiano (1988) to conclude that 
the nuclear emission was softer and more absorbed than 
that from the entire galaxy. A subsequent observation by 
Ginga in the 2-10 keV band (Ohashi et al. 1990) indicated 
that the spectrum was best fitted by a 6-7 keV thermal 
bremsstrahlung and no significant absorption; no iron line 
emission was detected (EW <400 eV). A more recent 
ASCA measurement (Ptak et al. 1997) required a two- 
component model to explain the complexity of the NGC 
253 X-ray spectrum [i.e., a soft thermal component with 
kT ~ 0.8 keV, plus a hard absorbed (A H ~ 10 22 cm~ 2 ) 
component, described either as thermal (kT ~ 7 keV) or 
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power-law (r ~ 2)]: again, no iron line emission was de- 
tected (EW <180 eV). At higher energies, a 4.4 a detec- 
tion by OSSE has been claimed by Bhattacharya et al. 
1994 with a 50-165 keV flux of 3 x KT 11 ergs cm" 2 s"\ 
higher than that extrapolated from the Ginga/ ASCA ob- 
servation. 

2. The BeppoSAX observation 

NGC 253 was observed on Nov.29-Dec.2, 1996 (see Ta- 
ble 1). The source is clearly extended (and elongated) up 
to ~ 10' in the BeppoSAX image in both the 0.1-2 keV 
and 2-10 keV band. Analysis of the resolved emission is 
deferred to a separate paper (but see preliminary results 
in Cappi et al. 1998): here we only present the analysis 
of the unresolved nuclear emission. In the central 4' re- 
gion (typical extraction region for point sources used for 
the LECS and MECS instruments), we obtained a LECS 
count rate of 3.93 x 10~ 2 cts s _1 and a MECS count rate 
of 9.23 x 10 -2 cts s _1 . Background spectra, extracted from 
the standard blank-sky files provided by the BeppoSAX 
Science Data Center (SDC), contributed less than about 
5% and 15% at 2 and 6 keV, respectively: similar results 
were also obtained using background spectra extracted di- 
rectly from the instruments' fields of view. LECS data 
above 4 keV were excluded because of remaining calibra- 
tion uncertainties. LECS and MECS data were reduced 
using the SAXDAS v. 1.3.0 and XANADU packages, while 
the XAS v.2.0 package was used for the PDS data. LECS 
and MECS spectra were rebinned in order to obtain at 
least 20 counts per energy channel, to allow use of the x 2 
statistics. 



Table 1: Exposure times. 



Instrument 


En. range 


Obs. time 




keV 


sec 


LECS 


0.1-4 


54689 


MECS 


1.3-10 


113403 


PDS 


13-60 


51557 



No short or long term variability is detected from the 
present data in any energy band: a fact that can be sim- 
ply interpreted as evidence for the large size of the emit- 
ting region. The 0.5-2 keV and 2.0-10 keV fluxes ob- 
served by BeppoSAX are ~2.53 xlO -12 erg s _1 cm~ 2 and 
4.9 xl0 _12 ergs _1 cm~ 2 . These are roughly consistent 
with the observed ROSAT PSPC and ASCA fluxes (Ser- 
lemitsos et al. 1996; Ptak et al. 1997) if one takes into ac- 
count the different sizes of the source regions. The global, 
absorption-corrected, 0.5-2 keV plus 2.0-10 keV luminosi- 
ties are 7.72 x 10 39 ergs s _1 and 1.57 x 10 40 ergs s _1 for 
the soft and the hard component (see Table 3), respec- 
tively, for an assumed distance of 3 Mpc. 



The BeppoSAX PDS count-rate between 13-60 keV 
is 0.07±0.02 cts s -1 which corresponds to a ^2.5 a level 
detection if systematics errors (~ 10%) are taken into ac- 
count: thus the 3<T upper-limit to the PDS 50-165 keV 
flux is 7.1 x 10~ 12 ergs cm~ 2 s _1 , much lower than the 
flux claimed to have been detected by OSSE. 

In the following, AT Hga i=1.28 x 10 20 cur 2 (Dickey & 
Lockman 1990) is assumed. 



Table 2: Bremsstrahlung plus emission lines model. 



Element ID 
for Kq, line 


Obs. Energy 
keV 


EW 
eV 


Fe L, Ne IX/X 


o.95i ;°o! 


109+tt 


Si XIV/XV 


l-9lt°o:°04 




S XV 


2 42+ 04 


79± 3 3 4 4 


Fe XXV 


6 69+ ' 07 


or,q+89 



Notes: The value of the relative normalizations Alecs/ Amecs 
is 0. 641^03- Intervals are at 90% confidence for one interesting 
parameter. 

* Line energy fixed to the ASCA value (Ptak et al. 1997). 
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Fig. 1. Thermal bremsstrahlung fit to the MECS and LECS 
data; the lines stand out clearly on the continuum. 



3. Spectral analysis 

We first used a thermal bremsstrahlung continuum plus 
emission lines model to parameterize the line energies and 
intensities. The LECS and MECS spectra were fitted si- 
multaneously, with free relative normalization. The best- 
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Fig. 2. Variable-abundances Mekal model fit to the MECS, 
LECS and PDS data. 

fit value for the central temperature is kT — 7A0t_Q'jl 
(with x 2 /dof=254/232). A few emission lines are evident 
in Fig.l: their resulting intensities and energies are listed 
in Table 2. In agreement with the ASCA results, most lines 
are consistent with K Q emission from H-like and He-like 
ions. 

At low energies 3 keV) the strong emission lines, 

notably the residual broad excess at E < 1 keV (a clear 
signature of the Fe L/Ne complex), require: (a) a soft ther- 
mal component with kT <1 keV; and (b) that the hard 
component be absorbed in order not to overproduce the 
continuum for Elz 0.7 keV (see also Ptak et al. 1997). 

At higher energies, an Fe K line at ~ 6.7 keV is un- 
ambiguously detected for the first time in NGC 253. It is 
narrow (a < 0.13 keV), with EW = 330 eV and consistent 
with emission from Fe XXV. [The measured equivalent 
width is roughly consistent with the upper limits placed 
by previous studies (Ohashi et al. 1990; Ptak et al. 1997) if 
one considers that those authors used a flatter power-law 
continuum.] Moreover, compared with ASCA, the larger 
MECS effective area at E>7 keV allows us to firmly estab- 
lish that the 2-10 keV continuum is dominated by thermal 
emission, since a hard power law alone (still allowed by the 
ASCA data: Ptak et al. 1997) yields a significantly worse 
fit (Ax 2 =20) with the present data. 

A two-temperature plasma is also required by the fol- 
lowing consideration. An atomic line energy is a function 
of the electron temperature kT and the degree of ioniza- 
tion (measured by n e t, where n e is the electron density 
and t the elapsed time). The centroids of the identified 
K-lines (Fe, Mg, Si, S: see Table 2)), imply allowable kT- 
n e t regions (see Fig. 2b of Kaneda et al. 1997) that are 
not compatible with a single-temperature plasma. Thus, 
at least two plasma components (with different ionization 
states) are required in order to reproduce the observed 
SAX spectrum of NGC 253. 



Therefore our best-fit model consists of a two- 
component thermal model. Best-fit parameters are re- 
ported in Table 3 for a two-temperature Raymond- 
Smith model and a two-temperature Mewe-Kaastra- 
Liedah plasma model ("mekal" in XSPEC). In order to 
reduce the number of free parameters, the abundances of 
He, C, and N were fixed at the solar value. The heavier 
elements were divided into two groups: Fe and Ni (most 
likely associated with SNe I products), and the a-elements 
O, Ne, Mg, Si, S, Ar and Ca (most likely associated to SNe 
II products); elements in the same group were constrained 
to have common abundances in solar units. 

Both components have consistent features within the 
statistical errors. The detection of the Fe K line in NGC 
253 allows us to reliably determine the Fe abundance in 
the line-emitting gas in the hard component: we find a 
value of ~ 0.25 solar, again consistent with the sub-solar 
values (based on upper limits) predicted by Ohashi et al. 
(1990) and Ptak et al. (1997). The iron abundance of the 
soft component is also constrained to a similar value. In- 
stead, the abundances of the a-elements are consistent 
with a single value of ~ 1.7 solar, for both spectral com- 
ponents. 

4. Discussion and Conclusions 

The likely emission mechanisms and environments for X- 
ray emission in SBGs have been described by Rcphaeli 
et al. (1995). In general, they may include thermal emis- 
sion from SNe, SNRs, and hot gas in the disk and halo, 
and nonthermal emission from massive X-ray binaries and 
Compton scattering of relativistic electrons off the FIR 
and CMB fields in the disk and the halo. Combined with 
the extent of the source of the 2-10 keV emission, the 
detection of the spectral lines reported here is a strong 
evidence for thermal hard X-ray emission from hot gas 
in SNRs and/or SN-driven wind, and severely limits any 
contribution from a nuclear low-luminosity AGN. 

A remarkable result from the present analysis is the 
striking similarity between the BeppoSAX spectrum of 
NGC 253 and the ASCA 0.6-10 keV spectrum of the 
Galactic Ridge X-Ray Emission (GRXE) in the Scutum 
arm region (Kaneda et al. 1997), i.e. unresolved diffuse X- 
ray emission components distributed in a thin disk along 
the Milky Way. The derived temperatures of both ther- 
mal continua and the identified lines with their centroid 
energies are similar for NGC 253 and the GRXE. 

Following Kaneda et al.'s detailed analysis of the 
GRXE soft component, the departure from ionization 
equilibrium that can be similarly inferred in the soft com- 
ponent of NGC 253 suggests that such component may be 
due to SNR-related emission from a diffuse plasma. For 
an estimated SNR luminosity of <~ 2 x 10 35 erg s _1 (see 
Kaneda et al. 1997), it can be estimated that the soft com- 
ponent's 0.5-10 keV luminosity is powered by ~ 4 x 10 4 
SNRs. The lack of appreciable photoelectric absorption, 
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implying a low HI column density where the soft com- 
ponent originates, also supports this galactic supcrwind 
picture. The derived chemical abundances of the soft com- 
ponent (see Table 2), that agree with the predictions of 
superwind and SBG chemical evolution models (Suchkov 
et al. 1994; Bradamante et al. 1998), favour too the galac- 
tic wind picture. In such models the galactic wind is pow- 
ered mostly by Type II SNe and minimally by Type la 
SNe: since Type II SNe eject mostly a-elements and Type 
la SNe eject mostly Fe, the resulting wind is overabun- 
dant in a-elements and underabundant in Fe. It should be 
noted, however, that uncertainties in the absorbing HI col- 
umn density and the likely spread of plasma temperatures 
may substantially affect the derived chemical abundances 
(see Read et al. 1997; Dahlem et al. 1998). 

The origin of the extended, thermal hard component 
is less clear. The substantial photoelectric absorption ob- 
tained from our best-fit requires a high in-situ hydrogen 
column density, Nh ~ 1-3 x 10 22 cm~ 2 (also the nuclear 
radio emission suggests the presence of molecular gas with 
similar column density: Nakai et al. 1987; Paglione et al. 
1996). This suggests that the hard component may orig- 
inate in buried, actively star-forming regions, with Type 
II SNe and X-ray binaries possibly supplying the bulk of 
the X-ray emission and the ambient gas and dust pro- 
viding the intrinsic absorption. In any case, this thermal 
plasma cannot be in equilibrium as it is not confined by 
galactic gravity. To estimate the contribution of SNe to 
the observed thermal flux, assuming that Type II SNe 
dominate the heating of the ISM, we repeat the swept-up 
mass calculation of Ptak et al. (1997) and similarly obtain 
that ^ 50% of the hard component originates in super- 
wind emission and not in the ambient ISM. The chemical 
abundances of the hard component, consistent with those 
deduced for the soft component, are likewise in agreement 
with the trend expected for SB-driven supcrwind models 
and hence are compatible with the SN origin for (a sub- 
stantial fraction of) the hard component. Moreover, Fe 
could be further depleted (relative to S and Si) by NGC 
253 's plentiful dust and warm ISM clouds (see Telesco 
1988). 

However, the galactic wind interpretation of NGC 
253's hard component raises some basic questions: (i) 
what is the mechanism leading to the high ~ 6 keV tem- 
perature?, and (ii) how many SNRs are required to power 
the hard emission? In fact: (1) SNRs usually have sub- 
stantially lower temperatures (fcT<4 keV); (2) even ne- 
glecting the spectral discrepancy, if ~ 2 x 10 35 
erg s _1 then 8 x 10 4 SNRs arc required to power 
the hard component of NGC 253, while usnr < 7.5 x 

10 4 (f^) 2 (o^)/!(fe) 3 ( where R SB and z SB are 
the radius and thickness of the SB region, and Rsnr is 
the SN shock front) SNRs can be at work if the two-phase 
status of NGC 253's ISM is to be preserved. Alternatively, 
the estimated ~4x 10 4 SNRs powering the soft compo- 
nent contribute only < 50% of the hard component; and 



(3) hydrodynamical simulations of superwinds do predict 
hard X-ray emission at such high temperatures, but with 
lower luminosities (compared to the soft X-ray compo- 
nent) and higher metal abundances (Suchkov et al. 1994): 
however, a self-consistent model of galaxy formation and 
evolution, including baryon infall and dissipation, star for- 
mation, and feed-back onto the still infalling gas, is re- 
quired to test the superwind picture (eg., Heckman et al. 
1990) in detail. 

Alternative explanations of the hard component in- 
volve either (1) a different mechanism leading to the ob- 
served thermal emission or (2) a "spectral conspiracy" 
(discussed in a separate paper) between a softer thermal 
emission and non-thermal emission yielding the observed 
spectrum. As for the former case, a magnetic confinement 
picture (e.g., Makishima 1994, 1995) entails that some 
fraction of the cooler plasma, confined by magnetic loops 
in the galactic disk, is heated to the observed temperature 
by the dissipation of the galaxy's rotational and velocity- 
dispersion energy through magnetic compression and re- 
conncction. As for non-thermal spectral components, the 
two most likely candidates are X-ray binaries and Comp- 
ton scattering of nuclear IR flux by relativistic electrons 
(Rcphacli et al. 1995; other nonthermal processes are dis- 
cussed by Valinia & Marshall 1998). As concerns point- 
like sources, Ptak et al. (1997), fitting ASCA SIS spectra 
for two strong extranuclear point sources by a Raymond- 
Smith plus power-law model, argue that the sources' fit 
parameters are similar to those for the entire galaxy, and 
the two points make up ~ 25% of the 2 — 10 keV flux. 
While some of the point sources in NGC 253 are associ- 
ated with SNe, some others are likely to be X-ray bina- 
ries with masses of at least 3 — 20 Mq (implied by their 
Eddington luminosities) suggesting that they are possible 
blackholc candidates. 

Summarizing, the main observational result reported 
in this Letter is the unambiguous detection in NGC 253 
of a hot luminous component whose nature is, however, 
still to be clarified. 
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Table 3: Two component-thermal models. 



Model kT so/t Ab so/t iv£ ard kT hard Ah hard X /do} 

keV a-el. xl0 22 cm~ 2 keV a-el. 

Fc, Ni Fe, Ni 

RS 0.88l:°; 22 1.52j:°;?° 1.28±g;|g 5.621°;^ < 0.80 246.6/237 

u - iy -0.05 u - zo -0.06 

Mekal 0.88l°;i^ 1.69±g;?l 1.30±g;tf 5.60t°;™ < 1.76 243.6/237 

O.lS^fr 0.27±°;J° 8 



Notes: The value of the relative normalizations A LE cs/Amecs is — 0.66 ± 0.03. Intervals are at 90% confidence for one 
interesting parameter. 
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